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We investigate the question of how tightly we can constrain the cosmological parameters 
by using the "cosmic inversion" method in which we directly reconstruct the power spectrum 
of primordial curvature perturbations, P(k), from the temperature and polarization spectra 
of the cosmic microwave background (CMB). In a previous paper, we suggested that it may be 
possible to constrain the cosmological parameters using the fact that the reconstructed P(k) 
does not depend on how many polarization data we incorporate in our inversion procedure 
if and only if the correct values of the cosmological parameters are used. The advantage of 
this approach is that we need no assumption regarding the functional form of P(k). In this 
paper, we estimate typical errors in the determination of the cosmological parameters when 
our method is applied to the PLANCK observation. We investigate constraints on h, Ob, 
Q m , and Qa through Monte Carlo simulations. 



§1. Introduction 

The cosmic microwave background (CMB) is one of the most powerful tools to 
determine the values of the cosmological parameters and examine the properties of 
the primordial fluctuations. A number of observations have been carried out since the 
Cosmic Background Explorer (COBE) observation. 1 ) In particular, a recent precise 
observation made using the Wilkinson Microwave Anisotropy Probe (WMAP) has 
confirmed that our universe is consistent with a spatially flat /1CDM universe with 
Gaussian, adiabatic, and nearly scale-invariant primordial fluctuations, as predicted 
by a simple slow-roll inflation model. 2 )~ 5 ) 

Nevertheless, it has been suggested that the primordial power spectrum of the 
curvature perturbation, P(k), may have some nontrivial features, such as a lack 
of power on large scales, running of the spectral index, and oscillatory behavior of 
the spectrum on intermediate scales. Therefore, a conventional parameter-fitting 
method, in which one often assumes a simple functional form of P(k), such as a 
power-law spectrum, is unsatisfactory. Instead we need to examine P{k) directly 
from observations without any theoretical prejudice. For this purpose, there have 
been several attempts to reconstruct the primordial spectrum using the WMAP data 
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using model-independent methods. 6 )~ 12 ) 

Cosmic inversion is a method to reconstruct the primordial spectrum directly 
from CMB anisotropics. It was originally proposed by Matsumiya et al. (2002, 
2003). 13 )' 14 ) We applied this method to the WMAP first -year data and showed that 
there are possible nontrivial features of P(k). 15 ^ Our method can reproduce fine 
features of P(k) with a resolution of Ak ~ 3.7 x 10~ 4 Mpc _1 , which roughly cor- 
responds to At ~ 5 in the angular power spectrum, Cg. In a previous work, we 
improved our method in such a manner that we can use the auto-correlations of 
both CMB temperature fluctuations (TT) and E-mode polarization (EE). 16 ) With 
polarization, we have shown that large numerical errors in the reconstructed P(k) 
due to the singularities in the inversion formula, which correspond to zero points of 
the transfer functions, can be suppressed. As a result, we were able to reconstruct 
P(k) with an error of a few percent in the ideal situation that observational errors 
do not exist. We have also found that it may be possible to constrain the cosmo- 
logical parameters by varying the contribution of the polarization in the inversion 
formula and requiring that the resultant P(k) is independent of the contribution of 
the polarization. In a conventional parameter-fitting method, it appears that one 
can determine the cosmological parameters up to a few percent from recent WMAP 
data. However, this is because the functional space of P(k) is restricted by the as- 
sumption of a simple functional form. As a result, these values of the cosmological 
parameters depend on this assumption. On the other hand, if we regard P{k) as 
a free function to be reconstructed, there remains the degeneracy that varying the 
shape of P(k) can compensate for the variation of the cosmological parameters. 17 )' 18 - 1 
Therefore, it is important to investigate how well we can determine the cosmological 
parameters without any assumption on the functional form of P(k). 

In this paper, we examine the proposition of determining the cosmological pa- 
rameters by using the cosmic inversion method. We first confirm that it is possible 
to determine the cosmological parameters quite accurately when there is no obser- 
vational error. Then we add artificial observational errors assuming the PLANCK 
observation*) and estimate probability distributions of the cosmological parameters 
by performing Monte Carlo simulations for each set of the cosmological parameters. 

This paper is organized as follows. In 321 we review our cosmic inversion method 
that employs both the CMB temperature and polarization spectra. We also extend 
it to a nonflat universe and examine the effect of observational errors on the re- 
constructed P(k). In 321 we describe a new method to constrain the cosmological 
parameters by using our cosmic inversion method and report the results of simula- 
tions to estimate the errors on the cosmological parameters. Finally, we present our 
conclusion in SJ3| 
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§2. Inversion method 



2.1. Formula 



First we present the formula to reconstruct P{k) using both the CMB tem- 
perature and polarization spectra. 13 )' 14 )> 16 ) For completeness, its derivation is de- 
scribed in Appendix^] We consider only scalar-type perturbations in which B-mode 
polarization is absent, and assume Gaussian and adiabatic primordial fluctuations 
throughout this paper. 

The angular power spectrum of the CMB anisotropy, C xx [where X and X are 
either the temperature fluctuations (T) or the E-mode polarization (E)], and the 
primordial power spectrum of the curvature perturbation, P(k), are related as 

c xx = l rdk k s mK xx (rjQ k)i (24) 
^ Jo K 

where K xx (rjQ, k) is the kernel specified by the Boltzmann equation. As we see, this 
is an integral equation for P(k). To solve it, we first tentatively adopt the following 
two approximations. One is the thin last scattering surface (LSS) approximation, in 
which we perform the time integration of the transfer functions within the thickness 
of the LSS. The other is the small angle approximation, in which we introduce 
the new variable r = 2c£sin(0/2), which represents the conformal distance between 
two points on the LSS. With these assumptions, we obtain a first-order differential 
equation for P(k) from the TT spectrum and algebraic equations for P{k) from the 
EE and TE spectra, respectively. However, they all have singularities corresponding 
to the zero points of the transfer functions that relate the primordial curvature 
perturbation to the temperature and polarization multiple moments. The presence 
of these singularities leads to large numerical errors near them in the reconstructed 
P(k), particularly when observational errors are taken into account. 

To avoid this difficulty, we construct a linear combination of the TT and EE 
formulas, introducing a free parameter a, as 

-k 2 f 2 {k)P'{k) + [-2k 2 f{k)f(k) + kg 2 (k) + akh 2 (k)] P(k) 

= S TT (k)+aS EE (k), (2-2) 

where f(k), g(k), and h(k) are time-integrated transfer functions within the thickness 
of the LSS, and the source functions are defined by 

OO -I o I tmax 



oo ^ 2i + 1 (I - 2)! Cf E ' obs / r 2 



s-w^jf*' E ^frf ^4-iH r ' (2 ' 4) 

— ^min & 

Here, C xx,ohs is the observed spectrum, and bf X '^ = cf X ' ex ^''/cf X ' a ' pv ^ > ' is the 
ratio of the exact spectrum to an approximated spectrum calculated from a fiducial 
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spectrum p(°\k), such as a scale-invariant one. This ratio b xx turns out to be 
almost independent of P(k), and it therefore plays the role of a corrector of the 
errors caused by the approximations. The boundary conditions are given by the 
values of P(k) at the zero points of f{k) as 

assuming that P'(k) is finite at k = k s . In Eq. (|2-2j) . the terms that have a prefactor 
of a come from the EE spectrum, and a controls the contribution of EE relative to 
TT. Because the positions of the singularities for TT and EE, given by f(k) = and 
h{k) = 0, respectively, are different, if we choose an appropriate value of a such that 
the contribution of EE is comparable to that of TT, the solution of Eq. (|2-2|) becomes 
numerically stable, even in the neighborhoods of the singularities. We found such 
an appropriate value of a to be in the range 10 13 — 10 15 if there is no observational 
error. 16 ) As shown in a previous work, because of the tight coupling of the photon 
and baryon fluids at the LSS, the transfer function for EE, h(k), is much smaller 
than those for TT, f(k) and g(k), by a factor of ~ 10~ 7 . The origin of the smallness 
of h(k) is briefly explained in Appendix [X] This is why the appropriate value of a 
is in the range 10 13 — 10 15 . 

Our original formalism described in Appendix [X] is based on a flat universe, 
where the spherical Bessel functions appear as the radial eigenfunctions of a Fourier 
expansion. Here, we extend our method to a nonflat universe. The dominant ef- 
fect due to the curvature of the 3-geometry is effectively absorbed by adjusting the 
conformal distance to the LSS so that the angular diameter distance in a curved 
background is properly recovered. This gives rise to a shift of the Doppler peaks, 
with the overall shape of C xx unchanged. The shape of C xx changes mainly on 
large scales (at small £), and its change on small scales is only a few percent. 19 )' 20 ) 
Thus, under the small angle approximation, we need only modify d in C« ' app de- 
pending on the curvature. This means that we can still use the spherical Bessel 
functions instead of the ultraspherical Bessel functions, which appear as the radial 
eigenfunctions in the curved background in cf X,apv . In fact, we find that using the 
spherical Bessel functions causes only a small error that can be corrected by b xx 
and does not affect the resultant P{k). 

In conclusion, for a nonflat universe, we need no modification of our inversion 
formalism described in AppendixlAl except for the adjustment of d, and the inversion 
formula (|2-2|) is applicable regardless of the geometry. 

2.2. Effect of observational errors 

Although in a previous paper 16 - 1 we showed that our method is effective if we as- 
sume that there is no observational error, it is also important to elucidate the effect 
of observational errors on the reconstructed P(k). The existence of observational 
errors causes numerical errors to be amplified near the singularities, as shown in our 
recent work, 15 ) where we used only the TT spectrum. Therefore, we reconstructed 
P(k) from C xx with observational errors, varying a, which represents the contri- 
bution of EE. We used the PLANCK observation and estimated the observational 
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errors, AC XX , by using the analytic formula 21 ) 

where / s k y is the sky coverage (/ s ky=l for a full-sky survey), iVpi x is the number of 
pixels, crpi x is the noise per pixel, o^eam is the beam size, and ' real \ s the real 
power spectrum. We set / s k y = 1, N p i x = 2.3 x 10 6 , <r p i x = 2.0 ^K and 3.7//K for 
TT and EE, respectively, and cr beam = 0.13° for the PLANCK 143GHz channel. 

First, for each £, we drew a random number from a Gaussian distribution with 
AC XX around a theoretical C xx , and then reconstructed P(k) from each simulated 
data set. Next, we estimated the mean and variance of the reconstructed P(k) at 
each k for 1000 realizations. Our numerical calculations for evaluating C xx and the 
transfer functions are based on CMBFAST*). We used Cf x up to £ m3jX — 1500 for 
both TT and EE. We assume a scale-invariant spectrum, k 3 P(k) = const., and set 
the cosmological parameters as h = 0.70, = 0.050, f2 m = 0.30, Qa = 0.70, and 
r = 0.20. In this case, the positions of the TT singularities are at kd ~ 70, 430, 680, 
1030, • • • , and those of the EE singularities are at kd ~ 230, 560, 860, 1180, 
where d ~ 1.36 x 10 4 Mpc. The range of the reconstructed P (k) is between the first 
and fourth TT singularities, i.e., 70 < kd < 1030, where kd roughly corresponds to 
I. To focus on the effect of observational errors, we assume that the cosmological 
parameters are known precisely. 

Figure □ shows the results for a = 0, 10 14 , 5 x 10 14 , and 10 15 . We see that the 
error in P(k) is very large near the TT singularities for a = 0, which is also the case in 
the situation studied in our recent work. 15 ) However, as a is increased, near the TT 
singularities the error is reduced, while near the EE singularities it is amplified. For 
a = 5 x 10 14 , the numerical error seems to be strongly suppressed near both the TT 
and EE singularities. On small scales (for I > 700), the overall error in P(k) becomes 
larger as a is increased, because the detector noise in Cf E becomes the dominant 
source of error for the PLANCK observation. In the absence of observational errors, 
we showed that P(k) is accurately reconstructed in the range 10 13 < a < 10 15 . 16 ) 
Taking the PLANCK observational errors into account, this range is narrowed due 
to the amplification of the numerical error, and the appropriate value of a is found 
to be ~ 5 x 10 14 , for which we can still suppress the numerical errors. 

§3. Constraining cosmological parameters 

3.1. Method 

As mentioned in a previous paper, 16 ) it is in principle possible to constrain the 
cosmological parameters in our reconstruction method as follows. We need to exam- 
ine the behavior of the reconstructed P(k) when the free parameter a introduced in 
Eq. I)2-2|) is varied. We find that as long as we use the correct values of the cosmolog- 
ical parameters, P(k) is accurately reconstructed for an appropriate value of a. On 
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Fig. 1. Reconstructed spectra with observational errors. We assume a scale-invariant P(k) and 
exactly known cosmological parameters, and we add the PLANCK observational errors to the 
theoretical Cj T and Cf E . The thin and thick curves represent the mean and variance, obtained 
from 1000 simulations, respectively. We plot them for four different values of a (0, 10 14 , 5 x 10 14 , 
and 10 15 ) which controls the contribution of EE relative to TT. The horizontal axis, kd, roughly 
corresponds to £. The TT and EE singularities are represented by the symbols □ and A, 
respectively. 

the other hand, if we use incorrect values, we obtain a particular deformed shape of 
the reconstructed P(k) that depends on the value of a. To be more precise, if we use 
a relatively large value of a, the deformation appears near the EE singularities, while 
if we use a smaller value, the deformation appears near the TT singularities. The 
point is that such deformations caused by the incorrect choice of the cosmological 
parameters indicate not only deviations from the actual P(k) but also inconsistent 
results for P(k) among different values of a. 

Let us denote the reconstructed P{k) for a certain value of a by P a (k) and 
introduce a quantity that represents the difference between P ai (k) and P a2 (k). We 
define 



where a± ^ a<i- From the above argument, we speculate that D takes its minimum 
value with respect to variation of the cosmological parameters at the correct values 
of these parameters. To confirm this speculation, we estimated the values of D for 
different values of the cosmological parameters, assuming that the values of these 
parameters used in H'2.2\ are the true values. We used Cf x up to £ max = 1000 
for both TT and EE to avoid the effect of a large observational error. Hence, we 
take /c m in and k max to be the first and third TT singularities, respectively. We used 
oi\ = 10 15 and 02 = 10 14 , because these values are appropriate to suppress large 
numerical errors due to the singularities. 

First, we assumed no observational error in C xx and estimated D as a function 
of the cosmological parameters. To see the dependence of D = D(h, f2b, ^m, ^a) on 




(3-1) 
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P(k), we performed calculations for three different shapes of P(k). We adopted a 
scale-invariant P(k), one with a large peak and dip, and one with a small oscillation. 
The spectrum with a peak and dip is expressed as 



k 3 P(k) = A i 1 + at exp 



(k - ktf 
o-i 2 



1 + a 2 exp 



(k-k 



2! 



\2 



°2 2 



-1 



(3-2) 



and we set a x = a 2 = 1, fci = 0.03 Mpc~\ fc 2 = 0.06 Mpc -1 , u\ = 0.01 Mpc -1 , and 
o"2 = 0.005 Mpc -1 . The spectrum with an oscillation is expressed as 

/ k 

1 + ao sin 



(3-3) 



and we set a = 0.1 and fco = 5 x 10~ 4 Mpc -1 . To elucidate the dependence of D 
on each cosmological parameter, we varied h, fi^h 2 , f2 m h 2 , and Qk = 1 — — ^yl 
one at a time, keeping the others fixed to the assumed real values. Explicitly, we 
varied h from 0.60 to 0.80 with Ah = 0.01, keeping Q^h 2 , f2 m h 2 , and fix fixed; Slf, 
from 0.040 to 0.060 with AQf, = 0.001, keeping h, Q m h 2 , and Qk fixed; fi m from 
0.20 to 0.40 with Af2 m = 0.01, keeping h, Q h h 2 , and i?^ fixed; and i?^ from 0.60 
to 0.80 with AQj\ = 0.01, keeping h, i??,/i 2 , and f2 m h 2 fixed. Thus, the number of 
grid points is 21 in each case. Note that we fix the optical depth r in our analysis, 
because it changes the shape of the spectrum only on large scales. 

From the results shown in Fig. [2 we find that regardless of the shape of P(k), D 
as a function of each cosmological parameter takes its minimum value at the correct 
value of that parameter in any case. With regard to the difference, the spectrum 
that has a large peak and dip tends to give larger values of D, while that with a 
small oscillation gives almost the same values as the scale-invariant spectrum. We 
also find that D is most sensitive to Qa or Qk- This is because the curvature affects 
the angular scale. It follows that the positions of the singularities, which correspond 
to the zero points of the transfer functions, shift to incorrect positions, and this 
causes large deviations of the reconstructed P(k). As described in this section, we 
have confirmed that our method to constrain the cosmological parameters is effective 
as long as the observational error can be ignored. This method is quite intriguing, 
because it requires no assumption regarding the functional form of P(k). 

3.2. Error estimation 

We also performed simulations to estimate the errors on the cosmological pa- 
rameters in the case that we observe CMB anisotropics by the PLANCK satellite. 
We generated 1000 realizations with the PLANCK observational errors and recon- 
structed P(k) from each realization in the same way as described in H2.2I For each 
realization, we calculated the value of D by varying the cosmological parameters 
and finding the minimum of D, which represents the location of the real values of 
the cosmological parameters, as discussed in the previous subsection. In practice, 
however, the location of the real values may be different from the minimum of D, 
due to observational errors. Therefore, we constructed histograms of the values of 
the cosmological parameters at the minimum of D from the 1000 realizations and es- 
timated their probability distributions by Gaussian smoothing. The assumed model 
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Fig. 2. D defined by Eq. I|3-l|l as a function of each cosmological parameter. In the separate plots 
in each panel, we assume P(k) as a scale-invariant spectrum, a spectrum with a peak and a dip, 
and a spectrum with a small oscillation. The true cosmological parameters are assumed to be 
h = 0.70, Q b = 0.050, n m = 0.30, Q A = 0.70, and r = 0.20. We vary h, tt b h 2 , fl m h 2 , and Q K , 
keeping the others fixed at the assumed true values, in the top left, top right, bottom left, and 
bottom right panels, respectively. (See the text for details.) 



is a scale-invariant P(k) with the same values of the cosmological parameters as used 
in H'2.21 Ultimately, we should perform a wide range parameter search in multiple 
dimensions. However, here we used limited ranges for the possible values of the pa- 
rameters and performed parameter searches only in one and two dimensions, mainly 
because the purpose of this paper is to carry out preliminary examination of how 
our basic strategy of unrestricting the functional form of P(k) affects the parameter 
estimation, and, practically, because our computations are quite time consuming. 

First, we varied each cosmological parameter individually, keeping the others 
fixed, as described in the previous subsection. The results are shown in Fig. |21 We 
find that the probability distributions are nearly Gaussian, and their peaks lie near 
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the correct values in all cases, as expected, since we drew random numbers from 
Gaussian distributions around the theoretical C? in our simulations, as mentioned 
in H2.21 For each cosmological parameter, we calculated its most probable value 
and the la error from its estimated probability distribution. The result is given in 
Table [I] We see that the most tightly constrained parameter is Qa or Qk , whose 
relative error is a few percent, provided that the other parameters are known and 
fixed. This is because D is most sensitive to the curvature, as mentioned in the 
previous subsection. 

To see possible degeneracies among the cosmological parameters, we also per- 
formed two sets of two-dimensional analyses. In these analyses, we varied Q\, and 
Q m , while keeping h and Qk fixed, and varied h and Qa, while keeping Q\,h? and 
Q m h 2 fixed. To save computational time, we investigated the same range of the 
parameter space, but with bin sizes twice as large as the corresponding ones in the 
one-dimensional analyses. Thus, the number of grid points is 11 x 11. The results 
for i?b and Q m are shown in Fig. 0J and the estimated values of the cosmological 
parameters are listed in Table [nj We find that in both cases the peaks deviate 
slightly from their correct values. In particular, the deviation from the correct value 
of fib is quite large. As shown in Fig. |SJ for h and Qa, there is a degeneracy caused 
by the fact that the same angular diameter distance can be the same for different 
sets of values of h and Qa 22 ^ We also see a peak near the correct values, but this 
peak is artificial, due to the sparseness of the parameter values we used. Moreover, 
the peaks in the projected probability distributions are amplified due to the narrow- 
ness of the parameter range, in which the global probability distribution does not 
converge sufficiently. For these reasons, we cannot estimate the values of the cos- 
mological paramters from this analysis. We need more fine-meshed and wide-ranged 
parameter search, although the computational time becomes quite long in this case. 

§4. Conclusion 

In a previous work, 16 ) we proposed a method to reconstruct the primordial 
spectrum by using both the CMB temperature (TT) and the polarization (EE) 
spectra, and we showed that it is effective if there is no observational error. We also 
proposed a new method to constrain the cosmological parameters by using the fact 
that the shape of P{k) obtained in this way must be independent of the contribution 
of EE relative to TT, which is controlled by the dimensionless parameter a. We 
found that the resulting P(k) depends very strongly on a, unless the correct values 
of the cosmological parameters are used in the reconstruction. Using this fact, we 
can, in principle, constrain the cosmological parameters without any assumption on 
the functional form of P(k). 

In this paper, first, to elucidate the effect of observational error, we have recon- 
structed P (k) from Ci with the errors expected from PLANCK satellite observations, 
assuming that the cosmological parameters are known. As mentioned above, the con- 
tribution of EE relative to TT is parameterized by a. We have found that numerical 
errors due to the singularities corresponding to the zero points of the transfer func- 
tions are suppressed for a ~ 5 x 10 14 , even when there exist observational errors. As 
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Fig. 3. Probability distributions of the cosmological parameters estimated with the one-dimensional 
analyses. The histograms were obtained from 1000 simulations, and the dashed curves are the 
estimated probability distribution functions. The assumed model is a scale-invariant spectrum 
with the same cosmological parameters as in Fig. H The bin sizes for the parameters are 
Ah = 0.01, An b = 0.001, Afl m = 0.01, and AQ A = 0.01. We vary h, Q b h 2 , tt m h 2 , and Q K 
individually, keeping the others fixed, in the top left, top right, bottom left, and bottom right 
panels, respectively. 

mentioned in a previous paper, 16 ) this value arises from the difference between the 
relative magnitudes of the transfer functions of TT and EE due to the tight coupling. 

We have also investigated the possibility of constraining the cosmological param- 
eters by using our inversion method, introducing the quantity D, which represents 
the difference between the reconstructed spectra for different values of a. In the ideal 
case with no observational error, we have shown that D takes its minimum value 
at the correct values of the cosmological parameters, independently of the shape of 
P(k). 

Then, to determine how tightly we can constrain the cosmological parameters in 
a realistic situation, we have performed simulations by taking the PLANCK observa- 
tional errors into account. We generated 1000 realizations and determined the values 
of the cosmological parameters which minimize the value of D for each realization. 
By constructing histograms of these values, we estimated their probability distribu- 
tions and la errors. However, to save the computational time, we have performed 
only one- and two-dimensional analyses. 

In the one-dimensional analysis, where we vary only one of the cosmological 
parameters at a time, with the others fixed at the assumed values, we found that 
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Fig. 4. Probability distributions of the cosmological parameters estimated with the two-dimensional 
analysis in (Qb , O m ) space. The assumed model is the scale-invariant spectrum with the same 
cosmological parameters as in Fig.|!2] The bin sizes for the parameters are taken as AQb = 0.002 
and AQ m = 0.02. The values of h and Qk are fixed at the assumed values. The top panel 
shows the two-dimensional probability distribution, where the solid, dashed, and dash-dotted 
curves represent the la, 2a, and 3er regions, respectively, and the bottom left and right panels 
show the projected probability distributions for Qi and Q m , respectively. 



their probability distributions are nearly Gaussian, with the mean values close to the 
correct values. We also found that Qa is constrained most severely, with variances 
of a few percent relative to the mean values. This is because the variation of Qa is 
equivalent to the variation of Qk for fixed Q m , which significantly affects the angular 
diameter distance to the LSS, and an incorrect choice of Qa displaces the locations 
of the singularities in such a way that the consistency of the shapes of the TT and 
EE spectra is lost. This leads to a large value of D. 

In the two-dimensional analysis, we have investigated two cases, one in which 
only Qf, and Q m are varied, with h and Qk fixed, and one in which only h and Qa 
are varied, with Q^h 2 and Q m h 2 fixed. In the latter case, there is a degeneracy in 
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Fig. 5. Probability distributions of the cosmological parameters estimated with the two-dimensional 
analysis in (h , Oa) space. The assumed model is the scale-invariant spectrum with the same 
cosmological parameters as in Fig. The bin sizes for the parameters are taken as Ah — 0.02 
and AQa = 0.02. We vary h and Qa, keeping Qbh 2 and Q m h 2 fixed. The top panel shows 
the two-dimensional probability distribution, where the solid, dashed, and dash-dotted curves 
represent the la, 2a, and 3a regions, respectively, and the bottom left and right panels show 
the projected probability distributions for h and Oa, respectively. The short-dashed curve in 
the top panel represents the degeneracy curve, where the angular diameter distance is the same 
as that of the assumed model. 



the (h , Qa) plane, because the angular diameter distance to the LSS, which deter- 
mines the scale of the spectrum, depends on both h and Qa- 22 ^ In our analysis, the 
sparse mesh and narrow parameter range caused an artificial peak in the probability 
distribution. Thus, we need to improve our computational method by optimizing 
the numerical code. 

In any case, from the one-dimensional analysis, we may conclude that even if 
we allow an arbitrary functional form of P(k), the cosmological parameters can be 
constrained by our inversion method, though this conclusion must be regarded as ten- 



Constraining Cosmological Parameters by the Cosmic Inversion Method 13 

Table I. Estimated values of the cosmological parameters obtained from the probability distribu- 
tions shown in Fig. [3] 





Assumed value 


Estimated value 


Fixed parameters 


h 


0.70 


n 7n o+0.027 

u.^uo_ 025 


fibh 2 , flmh 2 , Ok, r 




0.050 


04Q1+ 0031 


h, S7 m h 2 , Q K , t 




0.30 


n 2 oo+0.020 


h, fl b h 2 , (Ik, t 




0.70 


701+ 011 
u -' ul -0.011 


h, S7 b h 2 , S7 m h 2 , t 



Table II. Estimated parameters obtained from the probability distributions shown in Fig. 2] Note 
that the range of Q b we investigated may not be large enough, as seen in Fig. ^] 





Assumed value 


Estimated value 


Fixed parameters 




0.050 


0465+ 0037 

U.U4DO_ 0036 


h, f} K , t 




0.30 


n OOO+0.030 
u - zsz -0 .028 


h, n K , t 



tative, because we have explored only a limited range of the full, multi-dimensional 
parameter space in this paper. 

It should be mentioned that the constraints obtained with our method are, of 
course, weaker than those obtained with the conventional parameter-fitting method 
if the primordial spectrum P(k) exhibits a simple power-law form, as assumed in 
the latter. As the observational accuracy increases, however, there is a good chance 
that we will observe nontrivial effects of fundamental quantum physics 23 -* -25 -* and/or 
non-slow-roll evolution of the inflation-driving scalar field. 26 )~ 29 ) These effects may 
impart complicated features on the primordial power spectrum that are beyond a 
simple power law. If this is indeed the case, we would not be able to obtain accurate 
values of the cosmological parameters using the conventional method, which restricts 
the spectral shape from the beginning. Our method would serve as a powerful tool 
in such a situation and could be more appropriate for the next generation of high- 
precision experiments, such as PLANCK, which are intended to provide data from 
which we can derive not only accurate values of the cosmological parameters but 
also a more precise shape of the primordial spectrum. 
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Appendix A 

Derivation of Inversion Formula 



Here we review our method of the reconstruction of P{k) from both the CMB 
temperature and polarization spectra for scalar-type perturbations. 13 ^' 14 )' We 
assume Gaussian and adiabatic primordial fluctuations. Although we restrict our 
discussions to a flat universe, it is easy to extend our formalism to a nonflat universe 
as mentioned in Sec. 12.11 

The angular power spectrum of the CMB anisotropy is expressed as 

1 f°°dk k?(Xj(r)o,k)Xt{Tio,k)) 



4vr ° l ~2^Jo T 27TT ' (A ' 1} 

where Xe(r],k) and Xi(rj,k) are either 6^(77, /c) and Ei(rj,k) which represent multi- 
pole moments of temperature fluctuations and E-mode polarization in Fourier space, 
respectively, k is the comoving wavenumber and ry is the conformal time with r/Q be- 
ing the present value. These are expressed in the integral form of the Boltzmann 
equations (£ > 2): 30 ) 



0e(vo,k) _ r° 
2£ + l 



+ V b V(ri)j'e(kAr,) + \n 2 V( V ) [^{kAr,) + Je (kA V )] j>, (A-2) 



3 /(* +2 )i r dr?il2V(r?) M^ ; (A . 3) 



2^+1 2y (£-2)1 Jo yu (kA V ) 2 

where U2 = {02 — v / 6-E'2)/10, Arj = r]o — rj, and the overdot denotes a derivative 
with respect to the conformal time. Here V b is the baryon fluid velocity, 'F and F are 
the Newton potential and the spatial curvature perturbation in the Newton gauge, 
respectively, 31 - 1 and 

V(r])=fe- T ^, t(t?) = r°fd7], (A-4) 



are the visibility function and the optical depth for Thomson scattering, respectively. 
In the limit that the thickness of the last scattering surface (LSS) is negligible, we 
have V(rj) ~ 5(rf— ry*) and e~ T ^ ~ 0(rj—rj^), where 77* is the recombination time when 
the visibility function is maximum. 32 - 1 To obtain a better approximation, we take 
into account the thickness of the LSS. This is required especially for the polarization, 
since the CMB polarization is mainly generated within the thickness of the LSS. The 
approximation is to neglect the oscillations of the spherical Bessel functions in the 
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integrals. Applying this approximation to Eqs. I)A-2|) and (jA-3|) . we have 



2£ + l 



( r^dn [(0 O + ¥)V(ti) + (?- £)e- T W] | k (kd) 

Wr7*start J 



?7*start 



where d = % ~~ ^* is the conformal distance from the present to the LSS and ?7*start 
and 7?* e nd are the times when the recombination starts and ends, respectively. We 
have also replaced V& by @\ and neglected the quadrupole term in Eq. (|A-2|) . by 
adopting the tight coupling approximation. 32 ) We define time-integrated transfer 
functions within the thickness of the LSS, f(k), g(k), and h(k) as 

dn (0 o +^)( V ,k)V(n) + (^ -<P)(rj,k)e~ T{ri) = f(k)$(Q, k), (A-7) 

^7* start 

v * cnd d v G^r,, k)V(v) = s(*0*(o, fe), (A-8) 

?7*start 

-llj'» vwsl(iw0 ' l) ' (A ' 9) 

Here we have separated terms of the transfer functions which are dependent only on 
the cosmological parameters, and the primordial curvature perturbation which leads 
the primordial power spectrum is defined as 

P(k) = (\<P(0,k)\ 2 ). (A-10) 

Let us compare the magnitudes of the transfer functions. At the recombination the 
quadrupole, dipole, and monopole are related as 

& i ~ -tt( o + V) ~ M©o + 9) ~ kd ( ^ ) (0 O + SO, (A-ll) 
'' /2 1 k -01-7^—^1 -A f^Of— K (A' 12 ) 



(&<i) 2 (fca!) 2 a*n e aT (fed) 2 re e or fed \r]o / \n e a"T 

where n e is the number density of free electrons, o~t is the cross section of the 
Thomson scattering, and the subscript * denotes the value at the recombination. 
Since rj^/rjo ~ 0.02 and the mean free time of the Thomson scattering is much shorter 
than the cosmic expansion time, H^/(n e ax) ~ 10~ 3 , 31 ) 0i ~ 10~ 2 (kd)(&o + \P) and 
n 2 /{kd) 2 ~ lO-^W)- 1 ©!. Thus, at kd ~ £ ~ O(10 2 ), /(Jfe) ~ ^(fc) ~ 10 7 /z(/c). 

Substituting Eqs. (|A-5|) and (|A-6|) into Eq. (|A- If) . we obtain the approximated 
TT, EE, and TE angular power spectra, 

l^ TT ' apP = ^ r d T h * m [fWhW+gmw]** (Ai3) 
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21 + 1 cf E „„ = 2^+1 jf* *iW W *JiKH)]" , (A-14) 



4tt 



2*+l r TE,app = 21+1 / (I + 2)! 
4tt < 2^ 2 \j (£-2)\ 

00 dk 



o 



/c fe 3 P(AO [f(k)j £ (kd) + fc(fc)j<(foO-(A-15) 



The angular correlation function of the CMB temperature fluctuations is defined as 

~ of + 1 

C TT (#) = (0(ni)0(n 2 )} = V ——Cf T P e {cos8), cos9 = h x -h 2 . (A-16) 



47T 



Similarly, we introduce the following quantities for the polarization: 



^ TE W-E^\/|^f CPPrfcosS), (A-18) 

which are not the conventional angular correlation functions but they turn out to be 
convenient for inversion. Here we introduce a new variable r instead of 6 defined as 

Q 

r = 2dsin-, (A- 19) 

which is the conformal distance between two points on the LSS. From now on, we 
focus only on small angular scales, corresponding to r <C d, which is valid where 
^ £ O(10). For the TT spectrum, substituting Eq. (IA~T5|) into Eq. (lA~16l) . we can 
derive formulas for the approximated angular correlation functions in terms of P(k). 
With the help of the Fourier sine formula, we obtain a first-order differential equation 
for P(k), 

-k 2 f\k)P'(k) + [-2k 2 f(k)f'(k) + kg 2 (k)} P{k) 

= 4vr/ dr-^-{r 3 C TT ' app (r)} sin kr = S TT (k). (A-20) 
Jo r or 

This is the basic equation for the inversion of the TT angular power spectrum to 
the primordial curvature perturbation spectrum. Although the above differential 
equation is singular at f(k) = since the transfer functions are oscillatory, we can 
find the values of P(k) at such points, say k = k s , as 

p(m= S(M for f(ks)=0, (A ' 2i) 

assuming that P'(k) is finite at k = k s . We can then solve Eq. ()A-20|) as a boundary 
value problem between the singularities. Similarly, for the EE and TE spectrum, 
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substituting Eqs. (|A-14j) and t|A-15j) into Eqs. (|A-17j) and (|A-18|) . respectively, and 
so on, we obtain the algebraic equations for P(k), 

/•oo 

kh 2 {k)P{k) = 4vr / dr r C EE ' app (r ) sin kr = S EE (k), (A-22) 
Jo 

/•oo 

kf(k)h(k)P{k) =4x I dr rC TE ' app (r) sinAr = S TE (k). (A-23) 
•/ o 

In this case, we can find P(k) except for the singularities of h(k) = for EE, 
and those of f(k) = and h(k) = for TE, respectively. These are the basic 
inversion formulas for the EE and TE angular power spectra. If we use inversion 
formulas l|A-20[) . (|A-22|) and (|A-23j) separately, the reconstructed P{k) suffers from 
large numerical errors around the singularities of the respective equations. However, 
the fact that the zero points of f(k) and h(k) are different from each other can be 
used to resolve this numerical problem as follows. We construct a combined inversion 
formula by multiplying Eq. QA-22J1 by some factor a and adding it to Eq. (|A-20[) . 

-k 2 f 2 (k)P'(k) + [-2k 2 f(k)f(k) + kg 2 (k) + akh 2 (k)] P(k) 

= S TT (k)+aS EE (k). (A-24) 

The boundary conditions are similarly given by the values of P(k) at zero points of 
f(k) as 

"<■"■-> = Q^Sm f ° r «*■) =°> < A - 25 » 

assuming that P'(k) is finite at k = k s . Here we have taken a to be independent of 
k for simplicity and this free parameter a controls the contribution of EE relative to 
TT. If we take an appropriate value of a so that the contribution of EE is comparable 
to that of TT, the solution of Eq. (|A-24j) becomes numerically stable even around the 
singularities. This is because the contribution of EE dominates near the singularities 
of TT given by f(k) = 0, and vice versa. We have found such an appropriate value of 
a is ~ 10 13 — 10 15 if we assume no observational error. 16 ) The origin of this number 
comes from the fact that the transfer function h(k) for EE is intrinsically smaller 
than the transfer functions f(k) and g(k) for TT by a factor ~ 10 -7 as explained 
above and their squares are contained in the left-hand side of Eq. (|A-24|) , Since the 
TE formula, Eq. (JA-23|) . which is singular not only at h(k) = but also at f(k) = 0, 
is difficult to handle, we do not use it so far. 

For either TT or EE, the approximated spectrum (jXX,app ^ = j> OI ex _ 
pressed as Eq. (|A-13|1 or (|A-14[) has relative errors as large as about 20 — 30% to the 
exact one Cf X,ex expressed as Eq. (|ATj) . Unless we correct such errors due to the 
approximation, it leads wrong P(k) when solving Eq. (JA-24j) which uses the observed 
spectrum cf X ' ohs . Therefore, we introduce the ratio, 
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This ratio is found to be almost independent of P(k). We explain the reason as 
follows. In both C« ,ex and ' app , the transfer functions including the spherical 
Bessel functions are rapidly oscillating functions compared to a reasonable P{k) 
for sufficiently large i. Thus, when we take the ratio, the dependence on P{k) 
effectively cancels out. We use this property of b xx with some fiducial primordial 

spectrum P(°)(fc) such scale-invariant one. That is, we first calculate b« = 
(jXX, ex(o) j(jXX, app(o) ^ f rom p(o) anc [ then we use C XX ' ohs /b XX ' which gives 

a better guess to Cf X,app , instead of C XX,ohs , in the right-hand side of Eq. (|A-24|) , 
As a result, we obtain P{k) with much better accuracy. 

Finally, substituting Eqs. (|A-16|) and (|A-17|) into Eq. (|A-24j) and replacing Cf X ' app 
by C XX,obs /bf X ' , we obtain the inversion formula (|2-2j) in Sec. 12.11 
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